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Abstract 

We analyse the effects on the predictions for the microlensing searches toward the 
Galactic bulge coming from the fact that not all the stars monitored belong to the 
bulge itself, but that a non-negligible fraction of them actually are in the Galactic 
disk. The different distribution and motions of these disk stars make their associated 
microlensing rates and event duration distributions to be quite different from those 
of the bulge stars. We discuss the uncertainties in these predictions associated to the 
modeling of the Galactic components and the main implications resulting from the 
inclusion of this second source population. 
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Three experiments, MACHO (Alcock et al. 1995), OGLE (Udalski et al. 1994) and 
DUO (Alard et al. 1995), are actively searching for microlensing events in the direction of 
the Galactic bulge, and their observations arc already providing important information on 
the morphology of the inner Galaxy. Essentially two lensing populations are believed to be 
responsible for the events: the objects in the disk and those in the bulge itself (Griest et 
al. 1991, Paczyhski 1991, Giudice et al. 1994, Kiraga & Paczynski 1994). The large rates 
observed are most probably the result of the barred shape of the bulge (Kiraga & Paczyhski 
1994, Paczyhski et al. 1994, Zhao et al. 1995, Han & Gould 1995), with the major bar axis 
making a small angle with respect to the line of sight to the Galactic centre. A disk close to 
maximal would also be helpful in explaining the large rates observed. 

Up to now, an assumption made in all theoretical analyses is that the source population 
consists of the bulge stars alone. However, it is known that, in the fields observed, a non- 
neghgible fraction of the source stars actually belong to the disk. For instance, Terndrup 
(1988) estimated that 15% of the red giant stars in Baade's Window, at 6 = —3.9° and 
/ = 1°, belong to the disk. Probably this fraction is even larger for main sequence stars in 
the same field, due to the younger age of the disk. The fraction of disk stars should increase 
in fields at larger longitudes and similar latitudes, where microlensing observations are also 
carried out. It is our purpose in this letter to discuss the possible implications of this second 
source population for the analyses of microlensing results. 

Let us assume that, in a given field, a fraction / of the monitored stars belong to the disk, 
while 1 — / is the fraction of bulge stars. We denote tis the optical depth to microlensing 
of disk or bulge sources (s = D, B) arising from disk or bulge lenses {I — D,B). Hence, the 
total optical depth is just 

r = f[rDD + tbzj] + (1 - /) [tdb + tbb] . (1) 

We recall that the optical depth is averaged over the source distribution (Kiraga & Paczyhski 
1994) 

Tis^— dD^s-j^ / dD^i^pMD,, - D,i), (2) 

JMs Jo dJJos Jo C^J^os 

where Dos and Doi are the observer distances to the source and the lens respectively, pi is the 

lens mass density, dUg/dDos oc psD^^^^ describes the number density profile of detectable 

sources along the line of sight (Kiraga & Paczyhski 1994), and the normalization factor 

is Ng — /o^"""^ dDosdUs/dDos- The parameter (3 arises because the fraction of sources with 

luminosities larger than L is assumed to scale as L"^ . A reasonable range has been estimated 

to be /9 = 1 ± 0.5 in Baade's Window (Zhao et al. 1995), valid for 4 kpc< Dos < 12 kpc. 

We will adopt in our discussion j3 = 1, using this value also for nearby sources in the disk 

{Dos < 4 kpc), for which actually /3 should be somewhat smaller, since we do not expect the 

main results to change significantly with this last assumption. 

We will adopt for the disk distribution a double exponential profile with constant scale 
height 325 pc (Bahcall 1986), scale length 3.5 kpc and local density, in lensing objects, 
Pd = O.IMq/pc^, assuming hereafter that the solar galactocentric distance is Ro = 8.5 kpc. 
We will show results for Dmax = 12 kpc and 6 kpc, this last value in order to illustrate 
the possible effects of having a 'hollow' disk, as could be suggested by a certain deficit 
of disk stars beyond Dgg — 3 kpc inferred from the color-magnitude diagram (CMD) in 
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Table 1: Contributions to the optical depth (x10~®), rates Vis (in events/10^ stars yr) 
and average event durations (T) (in days), for lens and source populations in the disk (D) or 
bar [B). For F and (T) we assume lenses of mass m — Q.2Mq and use the OGLE efficiency. 
The first three columns are for Baade's Window, the last three for I — +10° and |6| = 3°. 

observations toward the bulge (Paczyhski et al. 1994b). We note that these observations 
may also suggest that the scale height is smaller toward the centre (see also Kent et al. 
1991). Finally, measurements of the scale length provide values in the wide range from 1.8 
to 6 kpc (Kent et al. 1991), further increasing the uncertainties associated to disk models. 
However, the hollow and non-hollow disk models considered already span a representative 
range of possible expected results. 

For the bulge stars, we will use the bar model of Dwek et al. (1995), obtained from 

COBE-DIRBE maps of the Bulge. The existence of a bar was initially suggested by several 
bulge observations, including an asymmetry in the infrared surface brightness distribution 
and in the star luminosities at positive and negative longitudes, and also to explain the non- 
circular gas motions observed (Spergel 1992). We adopt a total bar mass Mjjar — 2x IO^^Mq, 
in the upper range of different dynamical estimates (see however Blum 1995), as suggested 
by microlensing observations. The optical depths and rates due to bar lenses are of course 
proportional to Mhar, while the event durations are insensitive to it. The angle between the 
bar major axis and the direction to the Galactic centre is taken to be a = 20°, in the middle 
of the range a — 20° ± 10° obtained by Dwek et al.. We take the same velocity dispersion 
of bar objects as in Han & Gould (1995). 

In the third column of Table 1 we give the predictions for the different components of 
the optical depth toward Baade's Window. Since the optical depth of disk sources, td = 
'Tdd + tbd, is always smaller than the one for bulge sources, tb = t^b + tbb, taking into 
account the disk sources will lead to a smaller theoretical prediction for the total depth 
than in the case f — 0. This in turn will imply an underestimate, as noted by Bennet 
et al. (1994), of the required r^, and hence of the inferred bar and/or disk total mass 
normalizations required. For instance, if we assume / = 20% in Baade's Window, r is a 
factor 0.92 smaller than tb if Dmax = 12 kpc is considered, while it is a factor ~ 1 — / = 0.8 
smaller if D^ax = 6 kpc, since td turns out to be very small in this case. We also see that 
To is very sensitive to D^^x, but tb has a milder dependence on it which arises only through 
tdb- 
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Figure 1: Differential rate distribution in Baade's Window, assuming a common lens mass 
m = O.2M0 and a disk extending up to Dmax = 12 kpc. The curves are labeled by Is, where 
I, s = D, B for objects in the disk or in the bar respectively. 

We turn now to discuss the rates and event duration distributions. Of course, the rates 
have a similar decomposition as the optical depth in eq. (1). In Fig. 1 we show the differential 
rate for the different components as a function of the event duration, for Baade's Window 
observations. We assume Dmax = 12 kpc and take for definiteness all lenses to have a 
common mass m = 0.2Mq. 

A crucial difference between the disk and bulge sources, besides their spatial distribution, 
is their motion. Since disk objects have small velocity dispersion, here taken to be a = 
30 km/sQ, and a similar global motion due to rotation as the observer's one, DD type events 
are expected to have particularly long durations when Dos < Ro- These long duration events 
in Tdd, centered around 40 -7- 80 d, can clearly be identified in Fig. 1. There is a second 
contribution to F/j/) at smaller times (~ 5-i-30 d), due to events where the source star lies at 
Dos > -Ro) so that it is moving in the opposite direction than most of the disk lenses, which 
are at Doi < Rq. We assumed the circular speed of the disk to be everywhere 220 km/s, 
and we expect small modifications of the results obtained due to the slower actual rotation 
in the very inner bulge. It is apparent from Fig. 1 that the rates associated to disk stars 
are relatively more important at longer event durations, since short duration events are 
overwhelmingly due to Tbb (recall that the contributions to F from disk and bulge sources 
are weighted by / and 1 — / respectively, and that the observational efficiencies are generally 
very small for T < 10 d). 

In order to illustrate the dependence of the results in the assumptions done for Fig. 1, 
we show in Figs. 2 and 3 some variations of the assumed models. Fig. 2 considers the effect 

^the local dispersion of disk stars is ~ 20 km/s, but is larger toward the Galactic centre (see e.g. Lewis 
& Freeman 1989) 
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Figure 2: As in Fig. 1, but for bar models with different lens masses, m = O.IMq (BBl) 
and m — O.SMq (BBS) and for Kent's axisymmetric model with m — O.SMq (KK3). We 
also show the DD distribution for m — O.SMq. 

of changing the lens mass as well as the bulge model. There is no reason to expect that 
the mass function of disk lenses will be similar to that of bulge lenses. As an example 
to illustrate the implications that this may have, we show the predictions for BB lensing 
for masses m = O.IMq (curve BBl) and O.SMq (curve BBS), and compare them to the 
prediction of DD with m — O.SMq (curve DDS). If both masses are 0.3Mq, we have the 
same situation as in Fig. 1 slightly shifted to larger times (and with reduced rates). However, 
if the bar lenses have a mass smaller than the disk lenses, the relative contribution to the 
differential rate coming from the disk sources is more important at large times than before. 
Similar conclusions clearly follow if the total bar mass is actually smaller than the assumed 
2 X 10^°Mq. In the same way, if the triaxiahty of the bulge is reduced, the contribution 
coming from the bulge sources would be smaller than before. For instance, we show the 
predictions {KKS curve) for bulge-bulge lensing adopting Kent's axisymmetric model of 
the bulge (Kent 1992), obtained from Spacelab data, and assuming m = O.SMq. This model 
has a more centrally concentrated mass distribution and larger transverse velocity motions 
than the bar model, leading to shorter duration events and smaller rates. Similar results are 
obtained (Giudice et al. 1994, De Rujula ct al 1995) with the spherically symmetric heavy 
spheroid model of Caldwell and Ostriker, in which the bulge is just the inner spheroid. 
Although these models fail to account for the observed optical depth to the bulge. Fig. 2 
shows that any deviation of the bar model in the direction of the axisymmetric ones would 
tend to enhance the DD contribution at large event durations. 

Fig. 3 instead shows the effect of considering a 'hollow' disk. We depict, for m — 0.2Mq 
as in Fig. 1, the DD and DB contributions for Dj^ax — 6 and 12 kpc. Clearly the BD 
contribution, not depicted, becomes neghgible for Dmax — 6 kpc. We see that the long 
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Figure 3: Predictions for hollow {D^ax = 6 kpc) and non-hollow {D^ax = 12 kpc) disk 
models for Baade's Window. 

duration tail (T > 30 d) of the DB rates is not greatly affected, since it is mainly due to 
lenses at Doi < 6 kpc, far from the rapidly moving bulge sources. Instead, the DD rates 
become suppressed, with the short duration events clearly disappearing. 

In Fig. 4 we show the differential rate predictions at |6| = 3° and / = ±10°, for m = 
0.2Mq and Dmax = 12 kpc as in Fig. 1. The DD contribution is not significantly changed 
with respect to the situation in Baade's Window, except that the distinction between short 
duration events and long duration ones is less clear, due to the different transverse motion of 
disk objects at larger longitudes. The BB rates turn out to be less important in these fields. 
The rates involving the bar change significantly in the two fields, due to the inclination of the 
bar, and make Tdb quite important at negative longitudes, for the assumed D^ax = 12 kpc 
value, while Tbd (not plotted) becomes quite small. At positive longitudes, on the other 
hand, T^b is much smaller, is sizeable, and the DD contribution to the rate is the largest 
one for long duration events. Since the fraction of disk stars also increases with longitude, 
it becomes even more important to take the effects of disk sources into consideration in the 
microlensing searches in these fields. 

The fact that Tbd has the opposite behaviour than F^s in fields at positive and negative 
longitudes has the effect of reducing the asymmetric signatures in the microlensing maps 
of the bulge when disk source stars are taken into account. For instance, we obtain, for 
D^ax = 12 kpc and \b\ = 3°, that r(/ = 10°)/r(/ = -10°) ~ 1.5 if / = 0, while this ratio 
is ^ 0.94 if we take / ~ 0.4 (0.6) for / = 10° (-10°)g. We note that the asymmetry of 

^We note that a rough estimate of the angular dependence of the fraction /, in clear windows, may be 
obtained as // (1 — /) = kNd/Nb, where k is a constant chosen such that / ~ 0.2 at Baade's Window, and 
Nd and Nb are the values of Ns in eq. (2) computed for disk and bulge sources respectively. This leads to 
/ = 0.4 (0.6) for |6| = 3° and / = +10° (-10°). 
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Figure 4: Differential rates in fields at |6| = 3° and / = ±10°, assuming m = 0.2Mq and 

Dfnax — 12 kpC. 

the individual contributions between fields at positive and negative longitudes depends 
sensitively on the bar inclination assumed (and should disappear for a = clearly), but due 
to the above mentioned cancellations, it should be harder to get information on a from the 
asymmetries in the microlensing maps. However, the size of r do depend on a. For instance, 
the bar-bar optical depth at Baade's Window takes the values TBsioi = 30°) = 0.97 x 10~^, 
while tbb(q; = 10°) = 1.74 x 10~^, and hence the large rates observed suggest that the 
inclination is small. 

Going back to Table 1, we show in the last five columns the total rates and average event 
durations, assuming m = 0.2Mq and the observational efficiency of the OGLE experiment, 
for Baade's Window and for a field at |6| = 3° and / = +10°, as well as the optical depth in 
this last field. The average durations do not change significantly among the two fields, and 
reflect the behaviour of the differential rates just discussed in Figs. 1-4. From the total rates 
(including the efficiency), we may estimate the fraction F of events which result from disk 
sources, which is just F = fV^i/T . For Baade's Window, taking / = 0.2 we obtain F ~ 0.13 
if Dmax = 12 kpc, while F = 0.01 if Dmax = 6 kpc. In the field at |6| = 3° and / = +10° we 
obtain instead, assuming for illustrative purposes / = 0.4, that F ~ 0.39 if D^ax = 12 kpc, 
while F — 0.08 if D^ax — 6 kpc. So, it is quite plausible that some of the events among 
the almost 100 events observed up to now are due to disk sources. This fact is interesting 
also because it has been suggested that there is actually a possible excess of long duration 
events with respect to the predictions made with = I^db + ^bb alone (Han & Gould 
1995b). The long duration events associated to disk sources may help in this respect, since 
the total dV/dT distribution has an enhanced tail at large T. Finally, the fact that the time 
distribution of events is different when the disk sources are taken into account also affects 
the determination of the mass functions of the disk and bulge lens populations. 
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In the longer term, it may become possible to get, with enlarged statistics, more informa- 
tion about the microlensing from the disk source population, considering the lensing effects 
of stars in certain regions of the CMD. In particular, the bluer end of the main sequence and 
the brightest red clump stars, should in their majority be foreground disk stars (Paczyhski 
et al. 1994b). In the same way, considering only the region of the CMD where the bulge 
red clump stars lie, it is possible to minimize the contamination of disk stars (Bennet et al. 
1994), though never eliminate it[|. Moreover, the study of line of sight velocity dispersions of 
the lensed sources, and eventually of their proper motions, should also help to deduce their 
identity (Terndrup et al. 1995). Stars in a bar should have particularly large line of sight 
dispersion while the dispersion of disk stars should be much smaller. 

As a summary, we have considered the effects resulting from the fact that a fraction of 
the sources observed in microlensing searches toward the bulge belong to the disk. Since 
their optical depth is typically smaller than the one of bulge sources, the required tb needs 
to be larger than the measured optical depth, affecting then the model parameters inferred. 
This effect is maximum if the disk is hollow (assuming the same /), but is probably the 
only effect resulting from these disk models. If the disk is not hollow, the rates due to disk 
stars are also relevant. They are relatively more important for long duration events and will 
then modify the shape of the event duration distribution. The relative contribution from 
disk and bulge sources also depends on the mass function distributions of the two lensing 
populations, on the overall normalization of their densities and on the details of the bulge 
model. The contribution from disk sources increases with increasing longitudes, and can even 
become comparable to the one from bulge sources in the extreme field considered (/ = 10°, 
6 = 3°). For the barred model of the bulge, the longitude dependence of the rates has the 
opposite behaviour for disk sources than for bulge ones, reducing the overall asymmetries of 
the microlensing maps of the bulge. 



■^It is interesting to note that this cut seems aheady to lead to a larger optical depth than the standard 
cut including all sources (Bennet et al. 1994). 
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